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Abstract. In the present work, the non-spherical dust grain characteristics of 
comet Levy 1990XX with different silicate grain compositions (both pyroxene 
and olivine) are studied using the T-matrix method. Considering amorphous 
pyroxene grain (Mg^Fei-^SiOg) with x (= 0.4, 0.5, 0.6, 0.7, 0.8, 0.95, 1) and 
amorphous olivine grain (Mg2yFe2-2!ySi04) with y {— 0.4, 0.5) (Dorschner et al. 
1995), the observed polarization data of comet Levyl990XX are analyzed using 
the T-matrix code at A = 0.485 /xm. The best-fit spherical volume equivalent 
values of effective radius (re//), effective variance (we//) and shape parameter 
(E) obtained for the observed polarization data of comet Levy 1990 XX at 
A = 0.485/xm for amorphous pyroxene grains (Mga;Fei_a;Si03) with x = 0.4, 
0.5, 0.6, 0.7, 0.8, 0.95 and 1 are (0.213, 0.0015, 0.486), (0.215, 0.0022, 0.486), 
(0.215, 0.0019, 0.486), (0.217, 0.0034, 0.486), (0.219, 0.0057, 0.486), (0.222, 
0.0083, 0.486) and (0.230, 0.0092, 0.486) respectively. However, the amorphous 
olivine grains show a bad fit to the observed polarization data. Considering 
crystalline olivine grain in comet Levy 1990XX, Das & Sen (2006) (our earlier 
work) already studied non-spherical grain characteristics of that comet using the 
T-matrix method and found crystalline olivine prolate grains show satisfactory 
results. It is found that the shape parameter (E = 0.486) coming out from 
the present analysis for pyroxene grains is same with that of crystalline olivine 
grains. Thus the analysis of polarimetric data of comet Levy 1990XX at A = 
0.485 /im shows the presence of both olivine and pyroxene grains in that comet. 
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1. Introduction 

Comets are the least processed and the pristine materials of early solar nebula. Our 
knowledge of comctary dust comes from polarimetric studies of comets, remote observa- 
tion of IR spectral features and the in situ measurement of comets (e.g., comet Halley and 
the 'Stardust mission'). The polarization measurement of the scattered radiation gives 
valuable information about the shape, structure and sizes of the dust particles. Many 
investigators (Bastien et al. 1986; Kikuchi et al. 1987, 1989; Lamy et al. 1987, Le Borgne 
et al. 1987; Mukai et al. 1987; Sen et al. 1991a, 1991b; Chernova et al. 1993; Joshi et 
al. 1997; Ganesh et al. 1998; Manset & Bastien 2000, Das et al. 2004 etc.) have studied 
linear and circular polarization measurements of several comets. These studies enriched 
the knowledge about the dust grain nature of comets. 

Before the in situ analysis of cometary dust was possible in 1986, its composition 
was inferred from meteor spectra and laboratory analysis of interplanetary dust particles 
(IDP) thought to be related to comets (Millman 1977, Rahe 1981). But, the in situ 
measurement of Halley, gave us the first direct evidence of grain mass distribution (Mazets 
et al. 1987). Lamy et al. (1987) analyzed the data for comet Halley from spacecrafts 
Vegal, Vega II and Giotto. The important information about the chemical composition 
of dust particles in comet Halley has been obtained from the dust impact mass analyzer 
PUMA 1 and 2 on Vega and PIA on Giotto spacecrafts (Kissel et al. 1986a, 1986b; 
Mazets et al. 1987; Lamy et al 1987). The in situ measurement of comet Halley indicated 
three classes of particles: (i) The lighter elements H, C, N and O indicative of organic 
composition of grains called 'GHON' particles , {ii) Garbonaceous chondrites of Type I 
(G I chondrites) and (Hi) Mg, Si and Fe, called silicates (Glark et al. 1986). 

The infrared measurement of comets has provided useful information on the physical 

nature of cometary dust grains. Spectral features at 10 fim wavelength allowed the 
identification of silicates in comet dust. Another silicate feature at 20 ijm also appears 
to be present in many comets. The wavelength and shapes of these feature provide 
important information for the identification of the mineral composition (Wooden et al. 
1997, 1999; Hanner 1999; Hayward et al. 2000; Harker et al. 2002; Moreno et al. 2003 
etc.). 

The dust particles released by comets, are believed to contribute to the population 
of interplanetary dust particles (IDP), often get collected at high altitudes of the Earth's 
atmosphere. Laboratory studies have shown that majority of the collected IDPs fall 
into one of the three spectral classes defined by their 10 iim feature profiles. These 
observed profiles indicate the presence of olivine, pyroxene and layer lattice silicates. 
This is in agreement with the results obtained from Vega and Giotto mass spectrometer 
observations of comet Halley (Lamy et al., 1987). Mg-rich silicate crystals are also found 
within IDPs and are detected through cometary spectra (Hanner et al., 1999; Wooden 
et al., 1999, 2000). If these Mg-rich crystals are solar nebula remnants, then they must 
have formed in an environment close to the Sun (1-5 AU). Thus, the Mg-rich crystals 
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represent grains processed at high temperatures in the inner nebula. Therefore, Mg-rich 
grains in comets probe physical conditions and processes in the early solar nebula (Harker 
et al. 2002). 

In the present work, the observed polarimetric data of comet Levy 1990XX at A = 
0.485 /iTO (Chcrnova ct al. 1993) have been studied with different silicate grain composi- 
tions (both pyroxene and olivine) using the T-matrix code. 

2. Spheroidal grain model and the T-matrix theory 

It is now accepted that cometary grains are not spherical and may be fluffy aggregates or 
porous, with irregular or spheroidal shapes (Greenberg & Hage 1990). The measurement 
of circular polarization of comet Hale-Bopp (Rosenbush et al. 1997) reveals that cometary 
grains must be composed of non-spherical particles. In order to study the irregular grain 
characteristics of comets, T-matrix theory (Waterman 1965) is widely used by several 
investigators. Using the T-matrix code, Kerola & Larson (2001) analyzed the polarization 
data of comet Hale-Bopp and found prolate grains to be more satisfactory in that comet. 
Recently Das & Sen (2006) (our earlier work), using the T-matrix code, discovered that 
the prolate crystalline olivine grains can well explain the observed polarization in a more 
satisfactory manner as compared to the other shapes in comet Levy 1990XX. 

The T-matrix method is a powerful exact technique used to study the irregular grain 
characteristics of comet. In the present work, calculation has been carried out for ran- 
domly oriented spheroids using Mishchcnko's (1998) single scattering T-matrix code, 
which is available at http://www.giss. nasa.gov/'^ crmim. The main feature of the T- 
matrix approach is that it reduces exactly to the Mie theory when the particle is a 
homogeneous or layered sphere composed of isotropic materials. In the present work, the 
power law size distribution has been used for the analysis of polarimetric data, where the 
minimum and maximum particle radius are automatically set fixed for each and every 
run merely by specifying the particle effective radius (reff) and effective variance {veff)- 

The linear polarization of the scattered light by dust particles depends upon {i) 

wavelength of incident light (A), {ii) the geometrical shape (E) and size (r) of the particle, 
{Hi) Scattering angle, 9 (= 180°— Phase angle), and {iv) the composition of dust particles 
in terms of complex values of refractive index, m {= n — ik). It is to be noted that E > 1 
for oblate spheroids, E < 1 for prolate spheroids and E = 1 for spheres. 

The dust particles released by comets, are believed to contribute to the population of 

interplanetary dust particles (IDP). It has been observed that IDPs contain silicates both 
amorphous and crystalline in structure (Bradley et al. 1999). The major silicate species 
are pyroxene (MgajFei-^SiOs, where < x < 1) and olivine {Mg2y Fe2-2y Si04, where 
< 2/ < 1). A value of x = j/ = 1 denotes the Mg-pure end member of the mineral and a 
value of a; = y = denotes the Fe-pure end member (Dorschner et al. 1995). It has been 
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observed that comet Levy 1990 XX displays a strong silicate feature with a distinct peak 
at 11.25 fim, attributed to crystalline olivine grains (Lynch et al. 1992). 



In the present work, indices of refraction (n, k) of amorphous silicate minerals (py- 
roxene and olivine) (Dorschncr et al. 1995) have been taken. Since the refractive index 
parameters are not available at 0.485/um, so the values are taken at 0.50 /xm. It is to be 
noted that a higher Mg-content corresponds to a lower value of 'A:' (lower attenuation) in 
the UV, visible and NIR; and the more Fe in the mineral, the stronger the absorptivity 
in the UV, visible and NIR. Therefore, grains with less Fe and more Mg will be poor 
absorbers of solar radiation and cooler than their Fe-containing counter parts (Barker et 
al. 2002). Crystals in Chondritic porous interplanetary dust particles (CPIDPs) (Bradley 
et al. 1996) are Mg-rich {x < 0.9), solid and submicron in size (Bradley et al. 1999), and 
of probable cometary origin. 



The polarimetric data of comet Levy 1990 XX has been taken from Chernova et 
al. (1993). Since the polarimetric data is only available at A = 0.485/im, our analysis is 
restricted to that wavelength. Taking the grain composition of amorphous pyroxene grain 
with X = 0.4, the best fit values of rgg, Veg and E are determined using the T-matrix 
code. The calculations arc extended for x = 0.5, 0.6, 0.7, 0.8, 0.95 and 1 respectively. 
Now taking the grain composition of amorphous olivine (with y = 0.4 and y = 0.5), the 
calculation is repeated. The result coming out from the above analysis is shown in Table 
1. Also the result obtained from our earlier work (Das & Sen 2006) for crystalline olivine 
grain is also shown. 



It can be seen from Table 1 that pyroxene glasses with different grain compositions 
show significantly good fit to the observed polarimetric data of comet Levy 1990XX. It is 
interesting to note that the aspect ratio (E) coming out from the present analysis is same 
for all the grain compositions. Also, the effective radius of the pyroxene grains increase 
with the increase of x. In Fig 1, the simulated polarization curves with x = 0.4, 0.5, 0.6, 
0.7, 0.8, 0.95 and 1 are drawn respectively on observed polarization data of comet Levy 
1990XX. The; expected negative polarization curves have been successfully generated for 
different values of x using the T-matrix code. It can be inferred that the analysis with 
amorphous olivine grain shows bad fit to the observed data and is shown in Fig 2. It is 
important to note that the expected negative polarization curves have not been noticed 
with amorphous olivine grains. In Fig 3, the polarization curves for amorphous pyroxene 
(with X = 0.95 and 1) and the curve obtained by Das & Sen (2006) for crystalline olivine 
grain are shown. It can be seen form Table 1 that crystalline olivine grains show a 
slightly better fit to the observed polarimetric data as compared to amorphous pyroxene 
grains. 
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Figure 1. The simulated polarization curves generated using the T- matrix code at A = 0.485 
firn for amorphous pyroxene grain (Mg^ Fei-x SiOs) with x = 0.4, 0.5, 0.6, 0.7, 0.8, 0.95 and 1 
respectively. 

3. Discussions 



The study of the cometary spectra in the thermal region (8 - 45 /um) has revealed the 
presence of amorphous carbon and amorphous silicates, as well as crystalline olivine 
(Campins & Ryan 1989, Lynch et al. 1992, Brucato et al. 1999, Wooden et al. 1999, 
Harker et al. 2002 etc.). These findings are consistent with the composition of cometary 
dust particles as detected by Giotto and Vega spacecrafts (Kissel ct al. 1986a). The 
silicate emission features are also noticed in the observations of comet Hale-Bopp from 
ISO (Infrared Space Observatory) (Crovisier et al. 2000). Mg-rich silicate crystals are 
also found within IDPs and are detected through cometary spectra (Hanner et al. 1999; 
Wooden et al. 1999, 2000). Thus it is clear that Mg-rich grains are more satisfactory in 
comets. 



In the present study, the non-spherical grain characteristics of comet Levy 1990XX 
are studied using the T-matrix code, considering different silicate grain (amorphous and 
crystalline both) compositions. It can be seen from Table - 1 that Mg-rich amorphous 
pyroxene grains show good fit to the observed polarimetric data compared to amorphous 
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Figure 2. The simulated polarization curves generated using the T-matrix code at A = 0.485 
/Lim for amorphous olivine grain (Mg2;, Fe2-2y Si04 



Table 1. The size distribution parameter {tcb, Vcb) and shape parameter (E) at A = 0.485/^m 
for pyroxene (amorphous) and olivine (amorphous and crystalline) grains in comet Levy 1990XX. 



Material 


Refractive; index 


X 


y 


E 


rcir 


Veil 


A mm 




II 


/,■ 








(/"") 








1.709 


0.0515 


0.4 




0.486 


0.213 


0.0015 


9.93 




1.688 


0.0448 


0.5 




0.486 


0.215 


0.0022 


9.28 


Amorphous Pyroxene 


1.G78 


0.0291 






0.486 


0.215 


0.0019 


(i.54 


(Mg,Fei_,Si03) 


1.639 


0.0057 


0.7 




0.486 


0.217 


0.0034 


5.84 




1.613 


0.00179 


0.8 




0.486 


0.219 


0.0057 


5.46 




1..kS9 


0.00029 


0.95 




0.48G 


0.222 


0.0083 


5.30 




1.577 


0.00002 


1 




0.486 


0.230 


0.0092 


5.35 


Amorphous Olivine 


1.823 


0.104 




0.4 


0.461 


0.210 


0.0004 


82.30 


(Mg2;;Fe2-2j,Si04) 


1.768 


0.101 




0.5 


0.474 


0.213 


0.0019 


50.69 


Crystalline Olivine 


1.63 


0.00003 






0.486 


0.218 


0.0036 


5.22 


(Source: Das & Sen 2006) 
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Figure 3. The best fit polarization curves for amorphous pyroxene with x = 0.95 and 1 at A = 
0.485 fj,m. The curve for crystalline olivine grain from Das & Sen (2006) is also drawn here. 

olivine grains. Das & Sen (2006) already found that the crystalline olivine grains also show 

a good fit to the observed data. Thus it can be inferred from the present analysis that 
the dust grains in comet Levy 1990 XX are mostly amorphous pyroxene and crystalline 
olivine. It is interesting to note that the shape parameter coming out from the present 
analysis is same for both pyroxene and crystalline olivine grains and is given by 0.486. It 
has also been observed that the effective radii of the pyroxene grains increase with the 
increase of x (Ref: Table - 1). 

4. Summary 

1. The best- fit spherical volume equivalent values of rg//, We// and E obtained for the 
observed polarization data of comet Levy 1990 XX, at A = 0.485/im for amorphous 
pyroxene grains (Mg^jFei-ajSiOa) with x = 0.4, 0.5, 0.6, 0.7, 0.8, 0.95 and 1 are 
(0.213, 0.0015, 0.486), (0.215, 0.0022, 0.486), (0.215, 0.0019, 0.486), (0.217, 0.0034, 
0.486), (0.219, 0.0057, 0.486), (0.222, 0.0083, 0.486) and (0.230, 0.0092, 0.486) 
respectively. 

2. The expected negative polarization values have been successfully generated for pro- 
late pyroxene grains using the T-matrix code. 
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3. The amorphous ohvine grains show a bad fit to the observed data of comet Levy 1990 
XX and the simulated polarization curves do not show any negative polarisation. 

4. It can be inferred from the present analysis that the dust grains in comet Levy 1990 
XX are mostly amorphous pyroxene and crystalline olivine with same aspect ratio 
(E = 0.486). 
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